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ON Abstract 



We consider a directed random walk making either or +1 moves and 
CO a Brownian bridge, independent of the walk, conditionned to arrive at 

point b on time T. The Hamiltonian is defined as the sum of the square 
ON of increments of the bridge between the moments of jump of the random 

walk and interpreted as an energy function over the bridge configuration; 
£3 the random walk acts as the random environment. The thermodynamic 

i limit of the specific free energy is shown to exist and to be self- averaging, 

1 rt 

i.e. it is equal to a trivial — explicitely computed — random variable. 
O An estimate of the asymptotic behaviour of the ground state energy is 

also obtained. 



1 Definition of the model and main results 

Disordered systems are thoroughly studied nowadays as condensed matter models in the 
presence of impurities. Interesting questions, both mathematically and physically, concern 
the sample dependence or independence of intensive quantities like the specific free energy 
in statistical mechanics of spin glasses or the integrated density of states in spectral theory 
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or the critical exponents governing the asymptotic behaviour. Closely related models are 
random walks in random environments; in this context, interesting questions concern 
the possible modifications of the asymptotic behaviour critical exponents or the unusual 
normalisation required for the limit theorems to hold [12]. 

More recently, the methods developped in the statistical mechanical study of disor- 
dered systems have been applied to other systems like long protein or RNA molecules. 
It is argued in [6] that a protein molecule is very much like a random walk with ran- 
dom charges attached at the vertices of the walk; these charges are interacting through 
local interactions mimicking Lennard- Jones or hydrogen-bond potentials. On the basis of 
replica trick heuristics, it is claimed in [6] that the protein undergoes a phase transition 
between an unfolded state and a folded biologically active state. 

The purpose of this paper is to rigorously study a model arising in protein confor- 
mation. In [9], Kantor and Kardar introduced a model defined on the product space, 
Wjv X Qn, of iV-steps one-dimensional nearest neighbours random walks 

W N = { r : {1,--.,7V}->Z with r i+1 - r t G {-1, 1} for 1 < i < N - 1}, 

equipped with the uniform probability, with the probability space (Qjv, J 7 , F) where a 
collection q = (qi)\<i<N of symmetric, independent, and identically distributed random 
variables are defined. For given (r, q) } the Hamiltonian used in [9] is given by 

H5 K (r, q )= J2 m^, ri (i) 

l<i<j<N 

and has the interpretation of conformational energy of a protein molecule with N units 
whose geometrical shape is given by r and whose charges over the constituting units are 
given by q\ these charges interact through an ultralocal interaction. In [9], the random 
variables q are symmetric Bernoulli variables with values in { — 1,1}. This model is quite 
degenerate due to the discrete character of the q variables. For that reason, in [4], the 
model having the same Hamiltonian as that of the formula (1) but with the variables 
q being independent Gaussian was studied. Both models are however quite intractable; 
therefore, in [5] a further simplification was proposed: namely, the variables q are in- 
dependent Gaussian, the Hamiltonian is again given by the same formula (1), but it is 
defined on the space T>n of directed nearest neighbours random walks 

V N = {r : {1,...,7V}->N, with r l+1 - n G {0, 1} for 1 < i < N} 

instead of being defined on Wjv. This latter simplification allowed the use of transfer 
matrix techniques which, combined with numerical experiments, gave the asymptotic (for 

large N) behaviour of the ground state energy, inf H KK (r,q), of this simplified model 

rev N 

N 

under the constraint that the total charge, Q = is conditioned to remain fixed to a 

8 = 1 

constant b. 

Our aim is to go beyond the ground state energy and to study the thermodynamics 
of this model. The main difficulty however is to satisfy the constraint Q = 6, implying a 
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conditioning on this event, known to be a delicate procedure [10]. To avoid such subtleties, 
we decided to work on the space {Q = b}. Moreover, from a fundamental point of 
view, proteins are soft macromolecular structures, being able to vary their geometry 
although the charges of the constituting units remain fixed by the chemical properties 
of the molecule. Therefore, it seems more reasonable to consider the charges q as the 
microscopic configuration and the form (induced by r) as the environment. Our starting 
point is the Hamiltonian (1) considered as a function on the space Qn in a random 
environment r £ T>n. It proved computationally easier to consider the continuous version 
of the previous model; it is expected however that both discrete and continuous models 
have the same asymptotic behaviour. We stop here speculative considerations about the 
physical intrerpretation of our model and define it mathematically. 

Let T be a fixed strictly positive integer and b an arbitrary real number. We denote 
by Cjq' 6 t j the space of continuous real functions X defined on [0, T] such that X(0) = 
and X(T) = b. We denote by [X t )te[o,T] the standard Brownian bridge process on the 
space C^ T y Let (Yk)k>o be a sequence of independent Bernoulli variables, of parameter 
p = f{Yi = 1} = 1 — f{Yi = 0}, also independent of the Brownian bridge and defined 
on a probability space 0. Define, on the probability space 0, the renewal process, i.e. 
the strictly increasing sequence of positive integers (Ui)i>o } with Uq = and recursively 



Uk = inf{z > Uk-i '■ Yi = 1}. By Nj = sup{z : Ui < T} we denote the last index i such 
that Ui < T. Thus 

= U < *7i < • • • < U Nt < T < U Nt+1 <■■■ 
We set U = (U\, ■ ■ ■ ,Un) and S = (Si, • • • Sjv) with 



The energy of our model is given by the Hamiltonian Hj : X C — > M. , defined by 



This Hamiltonian is the continuous analogue of the Hamiltonian (1) considered in [9]. 
The length of the protein is identified with the parameter T and the Gaussian charges 
are modelised through the Brownian bridge process. The constraint of fixed total charge 
b is automatically satisfied by all configurations in C^ b T y The directed polymer acts as a 
random environment for the Brownian bridge. 

We define now the thermodynamics of the model. Denote by the mean expected 
value, by E{-|7£} the mean expected value conditioned to the process or variable 1Z, and 
by Ek{-} the mean expected value with respect to the distribution of 1Z. With these 
notations, the partition function is 



The parameter T plays obviously the role of volume. The finite volume, specific, quenched 
free energy is defined, as usual, by 



X Vi - X Vi 






(3) 



1 



logZ T (/3) 



(4) 
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and the finite volume, specific, annealed free energy is given by 



fAP) 



1 



logE( NtU )Z T (P). 



(5) 



This model, while reminiscent of the system studied in [11], has however a much richer 
structure since the random environment, given by the renewal process is independent of 
the bridge. 

With the formulation we use, we impose the total charge to be fixed, avoiding thus 
the usual subtleties [10, 1] of conditionning with respect to the final point. Moreover, the 
total charge being left as a free parameters from the very beginning, several scaling laws 
can be tested. 

We study the full thermodynamics of this model and we prove the following 

Theorem 1.1 Denote by V{ = Ui — Ui-\ the increment of the renewal process. Let 
lim 7j7 = z G [0,oo[. For every inverse temperature f3 G for every asymptotic 

T — > oo 

behaviour of the total charge z G [0,oo[ ; the quenched specific free energy, fx, converges 
in distribution, for T — >oo, to a trivial random variable f^, where 



andp = PfTx = 1). 

Remark: The previous theorem is established for a finite limit z. It is not difficult to see 
however that this result remains true even in the case z = oo provided that the divergence 
of the ratio b 2 /T is not very violent (sublinear). 

Besides the thermodynamic behaviour we obtain also results on the scaling of the 
ground state energy. It is argued, in [9], that when the total charge b scales asymptotically 
like b ~ T x (where / ~ g means limx^oo j{pj = 1) then the minimum of the energy scales 

like mini? ~ T e ^ x \ where 0(x) is a critical exponent continuously depending on x. These 
results are based on heuristic arguments and numerical simulations. Here we prove the 
following 



Theorem 1.2 Let < p = P(Fi = 0) < 1, q = 1 - p, b G R, and T G N + . Then, we 



foo(P) 



p 



Elog(l + 2PV 1 ), 



2/3 



have 




Corollary 1.3 Under the same conditions as the previous theorem 
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The previous result suggests that for fixed p and q = 1 — p } if b ~ T x , the dominant 
behaviour of the minimum of the Hamiltonian will behave as E[mini7] ~ T 2x . For p and 

X 

q depending on T however, the precise asymptotic behaviour can be changed. 

The paper is organised as follows: the next section deals with technical explicit com- 
putations and intermediate results on finite volume systems. These results are used in 
section 3 to obtain the thermodynamic limit of the specific free energy. The zero temper- 
ature limit and ground state energy are studied in section 4. In a last section some open 
problems are presented. The appendices remind some well known results for the paper to 
be self contained. 



2 The conditional expectation of the Boltzmann fac- 
tor 



Let u = (iti, • • • , u n ) be a vector with < u\ < ■ ■ ■ < u n < T . Denote by E(-|n, u) the 
conditional expectation with respect to the event {N = n; U\ = iti, • • • ,U n = u n } and by 
Hz = E(S|n,u) and Qz = (E[S 4 -Eij — (H2,)i(l J '2,)j\ n i u\; hi = = Gov = 

(Cov u} : z, j = 1, • • • , n) the mean expected vector and the covariance matrix 

of S conditioned to the event {(iV, U) = (n, u)} respectively. We have 

E[exp(-/3# T )|n,u] = (2tt)— / 2 |det Q^ 1 ' 2 j exp[-{/3# T + //^Ql 1 ^- He)}] \d£\. 

Denote q = 1 — p and let (j> : N + X W 1 -> 1 be an arbitrary integrable function; then we 
have [7] 



T-1 T-l-n T-l-(n-i) T-1 

E (N , u 0) = J2p n <i T - 1 - n E ••• E ••• E m(N,u)\ 

n=0 ui=l u t =u t -i+l ti„=u„_i+l 

Consequently, 

T-1 T-l-n T-1 

EZ T (P)='£p n q T - 1 - n J2 ■■■ E E[exp(-/3T/ T )|n,u]. 

n=0 ui=l u n =u n — i +1 



n,u) . 



Lemma 2.1 Set V{ = U{ — for i = 1, • • • , n and S 4 - = X Ui — X Ut _ 1 for i = 1, • • • , n. 
Denote by v = (v\, • • • , v n ) and for Vi ^ 0, denote by v~ x = (uf 1 , • • • , v n — l). Let D(y) 
the diagonal matrix given by D(v)ij = ViSij and, for every a £ IR ; denote by a the constant 
matrix whose coefficients dij = a for all i and j. Then, the inverse of the covariance 
matrix Qz exists and is given by 

Qz^Div-^+fr-p^ l = D(v- 1 ) + (T-u n )- 1 l. (6) 
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Moreover, 



det Q^ 1 = II vj 1 ( l + ( T - un)' 1 u n ) = n v ? )T(T- u^)" 1 • (7) 



Proof: From the elementary properties of the Brownian bridge reminded in appendix A, 
we get 

(fiz)i = E[Zi\n,u] = —Vi 

and 

Cov [(Ei,Ej)\n,u] = ViSij - —v.Vj for i,j = 1, ...,n. 

Introducing the matrix S(v) whose coefficients are given by S(v)ij = ViV 17 we can write 
the previous results in the matrix form 

fi E = and Q s = D(v) - ^S(v). 

n 

Since J2 Vi = ^ < 1, we can use lemma (C.l) to assert that Q E is invertible and more 

8 = 1 

precisely, obtain the explicit expression of the inverse 

Qi 1 = D(v~ 1 ) + ( r "X>) l = D(v- 1 ) + (T-u n Y 1 l. (8) 
Now use of lemma (B.l) allows to explicitely compute the determinant by the formula 



det Q^ 1 = J] vj 1 (l + (T - un)' 1 u n ) = J] v f }T(T- v^)' 1 



i=i \j=i 



□ 



Lemma 2.2 On the event {N = n; U\ = «!,••■,[/„ = u n }, the Hamiltonian can be 
expressed by the formula 

H T (u, E) = E'E + b 2 - 2bE f l + E'lE = E\I + T)E - 26S*1 + b\ 

where for a £ IR we put a = (a, • • • , a) £ W 1 to denote the constant vector of W 1 whose 
components are equal to a. 

Proof: On the event {N = n; U\ =«!,■••,[/„ = ti n }, the energy is given by Ht(u } X) = 

n n 

^(X^-X^Y + ib-X^) 2 . Since X Un = ES,-, we find 

8=1 8=1 



2 



n 



8 = 1 V 8 = 1 / 8 = 1 8 = 1 \8 = 1 
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Using the compact vector notation and denoting by (•)* the transpose of (•) (matrix or 
vector), we get the desired formula. □ 

Let us examine now the quadratic form 

M(vJ) = PH T (u,g) + fisYQi'd- //s). 

Since Q^ 1 is symmetric, we have 

M(u,i) = \(\Qz + 2/3/ + 2/31)f - i\2bfi+Q^^) + /36 2 + ^QsVs- 
Put Q(/3) = Ql 1 + 2/3/ + 2/31; from (8), we get 

Q(f3) = D(v)- 1 + 2/3/ + ((T - u n )- x + 2/3) T. 
With the compact vector notation introduced so far, 

g(/3) = /)(r 1 + 2/3) + 7 (/3)T 

with 

7(/3) = (T — Un)' 1 + 2/3. (9) 

Assuming that /3 > 0, the matrix Q(f3) is positive definite as it is shown in appendix 
B. Therefore, there exists a matrix T(/3) that is the square root of Q(f3). The notation 
T(/3) is an abridged writing of the more complete one, introduced in appendix B; namely, 

r(/3) = r(r 1 + 2/3, 7 (/3)); 

we stick however on this simplified notation not to excessively burden the formulae. Now, 
it is possible to express in terms of a vector r(/3) and a real c(/3) such that for any 

M(u,i) = l - (Y(P)i- r(/3)y (r(/3)f- r(/3)) + c(/3). 

Both existence and explicit formulae for r(/3) and c(/3) are obtained by remarking that 
T(/3) is a symmetric matrix and hence, by identification, T(/3)r(/3) = 26/31 + Q^ 1 [iz which 
implies 

r(/3) = r(/3)- 1 (26/31 + g^VH) (10) 

and 

c(/3) = fib 2 + ^QhVh " ^(/3)V(/3). (11) 

Lemma 2.3 Assuming ft > and using the notation introduced so far, we have 

E(exv(-PH T )\N,U)=ij>(N,U), 
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where 



ip(n,u) = E[exp(— (3Ht) \n, u] 



TiT-u^- 1 



1+7(/3)(e(^ 1 + 2/3)- 1 



> exp(-c(/3)). 



Proof: Writing explicitely the integrations involved in the conditional expectation, we 
obtain 



E 



[exp(-/3H T )\n,u\ = (2tt)-" /2 |det Q^e'^ J exp ( 



Make the change of variables £' = r(/?)£ — r(/3); this introduces a Jacobian in the volume 
element, of the form \d£\ = |det r(/?)| _1 |e?£'|. Then perform explicitly the Gaussian 
integration over the Brownian bridge to get 

E[exp(-PH)\n,u\ = |det Qi^^det r(/3)|" 1 exp[-c(/3)]. 
Now make use the corollary B.2 for the determinant of the matrix T(/3), namely 



|detr(/3)|= \J[(vf + 2P)j [l+mYX^ + m-'j ■ 

Combination with the formula (7) for the determinant of the matrix Q~_ leads to the 
desired result. □ 

Next we need an explicit evaluation of the parameter c(/3). This evaluation is quite 
straitghtforward but rather cumbersome; several auxiliary quantities are introduced that 
allow to express this parameter. Let 



2b 2 f3 2 



M{P) 
MP) 

MP) 

Am 



l(P) + [Eiv- 1 + 2/3)-i 



2b 2 (3 



t 



2(T - u 



2 Up) + (xx 1 + 



i=l 



\i=l 
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Lemma 2.4 With the notations introduced above, and under the conditions of the pre- 
vious lemma, we have 



c{(3) = (3b 2 + - 



1 b 2 U r . 



2 T T - u r , 



A 1 (/3)-A 2 (/3)-A 3 (/3)-A 4 (/3). 



Proof: Start from the expression (11) for c(/3) = (3b + ^^eQe P'Z ~ 2^(^)*^(^)' First 
evaluate the term ^[J-zQz fJ-E] using the expression given in (8) for the inverse of the 
matrix we get 



1 b 2 u n 



2 T T — u Ti 



Compute now the term -r(/3)*r(/3). From (10), the symmetry of the matrices, and the 
fact that T({3)- 2 = <5(/3) _1 , we get 

Ir(/3)V(/3) = 2b 2 {3 2 l t Q({3)- 1 l + 2b(31 t Q((3)- 1 Q^ V~ + ^sQl* W) _1 Ql Vh (12) 

Recall that Q({3) = D(v~ x + 2/3) + ^({3)1 with 7(/3) given by (9). Since we restrict 
ourselves to (3 > we can apply Lemma C.l to get 

Q(P)- 1 = D(ti(P)) - h(P)S(ti(P)), 

where 

w({3) = (v- 1 + 2/3)- 1 

and 

HP)= (7(/3)" 1 + E^(/3)) • 
Then the first term of the right hand side of (12) is equal to 

2b 2 p\i t DWP))i-h(p)i t sw))?) = 2b 2 p 2 ((£w0))-h(p)(jrw0)) ) 



\l = l 



26 2 /? 2 7 (/3)+ EK :1 +2/^)" 1 



\i=l 



so that 

2b 2 f3 2 l t Q(f3)~ 1 l = Ai(/3). (13) 

For evaluating the second term of (12), we first compute Q{fi)~ 1 Q^ 1 ■ On that purpose, 
introduce the following notation: for a vector £ denote by L(£) the matrix 

( £ x ■■■ £ x \ 

£ 2 ••• £2 



L(£) 



V 



/ 
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Observe that D(£)l = L(£) and S(£)l = (^J2 L($). Denote wv = (wiVi, w n v n ). 
Then we get, 

QiPY'Ql 1 = {D{w{P)) - h{p)S{w{p))} {Div- 1 ) + (T — u^ 1 T} 

= Diwip^Div- 1 ) + (T — u^ 1 (l - h{p) (^t «.■(/?))) L{w{p)) 
-hifiSWfflDiiT 1 ). 

Since 1 - h(P) (E «,M) = l(P)- 1 (l(P)- 1 + ( tu t (P))) 1 we get that 

, ^-i-,-n , Z((^ 1 + 2/3)- 1 ) 



Qi/Sy'Qi 1 = Ddv- 1 +2/3)- 1 v- 1 ) + 



(T - « n ) 7 (/3) ( 7 (/3)-i + EC"," 1 + 2/3)" 1 ) 



8 = 1 



The second term of (12) now reads 

2b 2 8 -, 

2b^Q{py 1 Q^^ = —E^QW-iQ-^ 

^2 



-f l t D(w(p))D(v- 1 )v + v 

J 1 (r-u n ) 7 (/3)^ 7 (/3)-i+ EMP)) 



l t 5'(^(/3)) J D(u- 1 )u 



7 (^)-i + E Wi{ p) 

\ 8 = 1 / 

We have the following identities 

n 

l t D(w(P))D(^ 1 )v = l t D(w(P))l = '£w i (P), 



8 = 1 



l t L(w(P))v= 5>,-(/3) \l t v= P>>8(/3) 

\i=l ) \i=l J 

and 



Combining, we obtain, 



2b 2 f3 f n \ 2b 2 /3 ( n ^ M/ 3 ) lu »-.E ! "^) 



2b^Q{py 1 Q^^ = ^ ( £ MP) ) +^ ( E MP) n — - 

8 = 1 
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so that, finally, the second term of (12) reads 

2bfi t Q(p)- 1 Q^ 1 ^ = A 2 (P) + A 3 (f3). (14) 
It remains to evaluate the last term of (12); rewrite, 

and 



= 1* + (T - Un)- 1 uj* 



From this equation and its symmetric couterpart 
1 , , 1 



U 2 (T - u n y 2 {l'Z>(u?(/3))l - hifj^Siwifj))!) 

1 



V (T - u n )~ 2 |5>,-(/3) - h(P) (5>,-(/3) 
\b 2 (T - u n )- 2 L/3) + (±(v^ + 2/3)" 1 ) 



so that finally 



^IQe 1 W) _1 QeVe = MP)- (is) 

□ 

Remark: Since explicit formulae for the conditional expectation with respect to the event 

{(N,U) = (n,u)} 

are obtained, the conditional expectation E(-|iV, U) is also explicitly known. In fact, if for 
an integrable function F, 

E(F\n,u) = <p(n } u) 

then 

E(F\N, U) = <f>(N } U). 

So that, if we replace in the above formulae n with Nj and U{ or V{ with Ui and V{ 
respecitvely, we obtain the corresponding expressions for the conditional expectations 
E(-\N, U). 
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3 Thermodynamic limit of the specific free energy 



In this section the T — > oo limit ol the specific tree energy will be evaluated. Some notation 
and well known results on the renewal process are necessarry. Since Uq = and U n = 
inf{k > U n -\ : Yk = 1, the variables Y{ being independent and identically distributed, 
the renewal process (Ui) can be written as a sum of independent geometric variabless 

n 

U n = K', where V{ = Ui — Ui-\. For every integer x > 1, we have that P(Vi = x) = 

8 = 1 

pq x ~ 1 , where p = W(Yi = 1) = 1 — q. Therefore, fj, = EVJ- = K Of special interest is 
also the random variable Nt = sup {A; : Uk < T}. Obviously, lim Nt = oo. Notice also 

T oo 

that although Nt is not a stopping time with respect to the cr-algebra cr(Yi } • • • , Yt), the 
random variable Nt + 1 is, on the contrary, a stopping time. 



Lemma 3.1 Almost surely, we have, 



lim — — 



p. 



Proof: The proof is quite standard (see [2] for instance) and relies on the strong law of 
large numbers. □ 



Lemma 3.2 The finite volume specific free energy, fj, is a random variable reading 



— Eioga^o + ^iogT- — 



log(T - U> 



1 



log 



2/3T 

1 Un t 



1 + ((T - U Nt )~ x + 2/3) ^(K- 1 + 2/3)" 1 
A(/3) 



8 = 1 



2(3T(T-U Nt ) (3T ' 
where A(f3) is a random variable such that 

E(A(/3)\n, u) = + A 2 (/3) + A 3 (/3) + A 4 (/3) 

and Ai, . . . , A 4 are the expressions given in the previous section. 



v 

+ T 



Proof: Recalling the remark on conditional expectations given at the end of the previous 
section and lemmata (2.3) and (2.4), the present result follows immediately. □ 

Studying the thermodynamic limit of fj is equivalent in studying the T — > oo limit 
in the above formula. Notice however that besides random terms whose limit must be 
computed, there are also deterministic terms of the form b 2 /T where b is the total charge 
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of the protein, supposed to scale with the size T of the molecule like b ~ T a . It follows 
that 

b 2 ( 1 if a = 1/2, 
z = lim - = if a < 1/2, 
1 { oo if a > 1/2. 

We do not fix the exact behaviour here. We only assume that z £ [0, oo] in the sequel. 

Lemma 3.3 Let (X n ) n >i be a sequence of independent identicaly distributed random 
variables such that nP(|Xi| > n) = o(l) and {T n ) n >\ be a sequence of positive integer 
valued random variables satisfying — > c in probability, with < c < oo. Then, denoting 
by S n = EILi Xi, we have 

J™o ("^ ~~ ^(Xi^ttx^Kn})^ = in probability. 

Proof: This is a standard extension of the strong law of large numbers and can be found, 
for instance, as theorem 5.2.6 of [3]. □ 

Proof of theorem 1.1: It is enough to study the T — > oo limit of all individual terms 
appearing in the expression for fx given in lemma 3.2. 

Let Xj = log(l + 2f3Vj) < 0. Because ~P(Xj > T) vanishes superexponentially fast, 
it follows that TP(|Xi| > T) = oil) and we can apply lemma 3.3 to show that the first 
term converges in probability, 

— £log(l + 2^) P ™ b ^E[log(l + 2/^)]. 

The second term trivialy vanishes in the limit. As for the third term, remark that 
1 < T — Un < T } therefore it also vanishes in the limit. For the fourth term, remark that 

(T - Un)' 1 + 2/3 < 1 + 2/3 and £ (Vf 1 + 2/3)" 1 < E Vj = U Nt < T, hence this term 

gives also a vanishing limiting contribution. The fifht term yields, by definition, z. 

The sixth term involves the study of the random variable Ax = T — Un. Obviously, 
F(Ax > k) = F([/jv T+ fc — Un t = 0). Now, since Ax < T, this random variable cannot be 
a plain geometrical one. In the limit T — > oo, it can be shown however (see formula 6.5, 
page 194 of [7] for instance) that F(Aoo > k) = q k ~ 1 , for k = 1,2, . . .. Hence, the sixth 
term yields 

1 b 2 U Nt z . 

lim — -ttt — — = — — — m distribution. 

T^oo2f3TTT -U Nt 2/3AOO 



Because (V; -1 + 2/3)" 1 > (1 + 2/3)" 1 , it follows that 

(n t \ - 1 

^((V- 1 + 2/3)" 1 < (1 + 2/3)iV T - 1 



,i=l 
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and since Nj — > oo this inverse sum vanishes in the limit. Hence the seventh term con- 
verges, in distribution, to — -, 2 o ■ 

For the next term we use again the lemma 3.3 to show that 



2b 2 /3 1 

~ 7 Ft 



-\ \ prob 



-2pzE[(V- 1 + 2P)- 1 ]. 



. i=l 



For the following term we obtain, 

J2 + 2j g Ar) -i _ 1 Y^MVf 1 + 2/3)- 1 ] 

Again, use the observations that 2/3 + A? 1 > 2/3 + 1 and E£ t i(K _1 + 2 /?)" 1 > iV T (l + 2/3)" 1 
so that the big denominator 1 + [• • •] reduces to 1 in the limit T — > oo. Subsequently, this 
term yields: 

-2z [(1 + 2(3 A^)' 1 - pEftVr 1 + 2/3)" 1 ]] in distribution. 



The same arguments can be used to study the last term 



2f3T(T-U NT y 



2/3 + A T 1 + ^(Vr 1 + 2/3)" 1 



,i=l 



Obviously, again, the term (j2i^i(Vi 1 + 2/3) ^ vanishes in the limit. It remains a 
contribution of the form 

z z 
+ 



2/3AO, l + 2/3A c 



Adding together all the terms, the contributions of all but the first one conspire to 
disappear. □ 



4 Finite volume ground state energy 

We can now consider the behaviour of the ground state energy. First notice that a 
standard Laplace's argument guarantees that 

- lim i log [E(exp(-/3# T )|(iV, U))] = inf H T (U,X). 

Proof of the theorem 1.2: the proof proceeds by examining how the individual terms in the 
expression of Ft ((3) = — \ log Zt(J3) = T/x(/3), where /x(/3) is given in lemma 3.2, behave 
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in the limit ft — >oo. Namely, we examine the behaviour when ft — > oo of the individual 
terms for the expression below: 



F T (ft) = Tf T (ft) 



1 Nt 11 

-gio g ( 1+ 2 W ) + ¥ iogr--io g( r- f /,) 



1=1 



1 Un t _ A^ft) _ A 2 (ft) _ A 3 (ft) _ A 4 (ft) 
2ft(T-U Nr ) /3 /3 /3 p ■ 



+ b 2 



It is easy to see that the sole terms remaining in the limit ft — > oo are the fifth one 
giving trivialy a contribution equal to b 2 — and the seventh one — namely, 



P^oo ft 



lim 



1 + (2/3A T )"i + [E!! 2 ! (1 + (2/3K-)" 1 )" 1 
Adding the two non vanishing contributions, we get 

b 2 b 2 
lim F T (ft) = b 2 T = 

and the theorem follows. 

Proof of corollary 1.3: This is obtained by an elementary summation: 



i + iV 



E 



1 



1 



T-l 

E 



T 



Therefore, collecting all terms, 



k + v- (n+i) = ^(i-? a 



□ 



E inf H T (U,X) = E lim F T (ft) 



v xec. 



[0,T] 



[3 oo 



pT 



:i - ? J 



□ 



5 Conclusion and open problems 

It is established that the model studied here has well defined thermodynamics in the sense 
that the infinite volume limit of the quenched specific free energy exists; moreover, this 
limit is a trivial random variable — i.e. almost surely a constant. This result shows the 
sample independence of the thermodynamic properties of the model at every temperature; 
a feature common to all short range disordered systems. 
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The large scale behaviour of the average ground state energy is also very instructive. 
Suppose first that the parameters p and q are fixed — independent of T — and non trivial 
i.e. < p, q < 1. Then, 

E[MH T (U } X)] ~ - for large T 

and it has a finite value only when z < +00. Now, a finite and strictly positive z means 
that the total charge scales as T 1 / 2 for large T, while the number of visited sites from 
the directed random walks behaves like pT. On the other hand there is no obstruction 
in choosing p and q depending on T. To illustrate what it may happen, suppose that we 
choose p to scale asymptotically as a/T, with < a < +00. Then 

b 2 

E[MH T (U,X)] (1 -exp(-a)) 

x a 

and this expectation is finite provided that the total charge remains finite. In the same 
time, the total number of visited sites does not diverge. Of course various other asymptotic 
behaviours can be obtained by choosing different scaling laws for the parameters p and q. 
The result we have obtained allows an exhaustive study of all these scaling behaviours. 

What remains an open problem for the moment is the statistical mechanics study of 
the system. Namely, for every realisation U of the random environment, we can define 
a random finite volume Gibbs measure, /ix, on the space C^ b T ^ as a measure absolutely 
continuous with respect to the standard Gaussian measure, 7^, of the Brownian bridge. 
This Gibbs measure can be defined through its Radon-Nikodym derivative 

^L(U X) = ex P(-P H T(U, x )) 

An ' Z T (f3) 

where Zt((3) is the partition function that depends also on U . It is not yet clear whether 
a Dobrushin-Lanford-Ruelle construction is possible for this model. Neither is it clear 
whether there is a phase transition or not, should it be in the DLR or the weak sense. 
These fundamental questions for the statistical mechanics formulation of the model are 
presently under investigation. 

A second problem that is still under investigation concerns the thermodynamics of the 
model defined for genuine random walks and not directed random walks. This destroys 
the explicit renewal process structure of the model and makes the computations much 
more complicated but not untractable. 

Summarising, inspired from the protein folding statistical mechanics problems, we are 
confronted with a new class of systems, reminiscent of both spin glasses and random walks 
in random environments. 
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A Identities for the Brownian bridge 



Let (W s : s > 0) be a Brownian Motion starting from and introduce the process 
(Y s :s> 0) by 



Y t = (T-t) j 



1 dW„ 



T-s 
o 

Then, the Brownian bridge between X = and Xj = b is expressed by 

X t = b^ + Y t foTte[0,T]. 

It is immediate to verify the following identities that are reported here only for complete- 
ness 

E(Y t ) = 0, E(X t ) = ^, E(X t - X s ) = ±(t - s) 

and 

„ /T ,„ T-t bH 2 (T-t) 

ny t 2 ) = —t, E(xi) = — + '—^t. 

For s < t, we have also 



T ^ v- s)),E((X t -X s f)=(^-^ 



E((Y t - Y s f) = ^— — — ^(T -(t- s)), E((X t - X s f) =(^-±)(t-s) 2 + (t-s)<t, 



and 



y aT ((x t -X s )) = (t-s)-^(t-s) 2 . 



Finally, for Si < s 2 < t\ < t 2} we have 

(s 2 - si)(t 2 - h) 



E((Y t2 - Y tl )(Y S2 - Y S1 )) 



T 



E((X t2 - X h )(X S2 - Y ai )) =(^-^j (s 2 - Sl )(t 2 - U), 

and 

Cov ((X t2 - X tl ), (X S2 - X S1 )) = -i( S2 - Sl )(t 2 - h). 



B Determinant and positive definiteness for a class 
of matrices 

For £ = • • • ,£ n ) G IR n denote by D(l) the diagonal matrix, D(£)ij = Ifiij for i,j = 
1, • • • ,n. For a G IR denote by a the constant matrix whose elements are a 8J = a for 
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z, j = 1, • • • , n. From [13], page 92, we get the following property, 



n n 



det (D(£) + a) = IK+«EII 

Consequently, if £j ^ for every j, then 



j = l 8 = 1 j = l 



det (D(/) + a) = \J[£A (^l + aE^J- ( 16 ) 



/4 0\ /a ••• a 

'• r : 

\ o 4/ \ a ■ ■ ■ a , 



i 



(17) 



Remark that the matrix (D(£) + a) as well as all its minorants 
(D(£) + a)^ 

for k = 1, • • • n are symmetric and have the same form. 

Lemma B.l Assume that £{ > 0, for all i = 1, • • • , n ond a > 0. Then, 

1. all minorants satisfy det (D(£) + a)( fc ) > /or A; = 1, • • • n, 

2. the matrix (D(£) + a) «'s positive definite, 

3. there exists a positive definite symmetric matrix T(£, a) that satisfies 

{D{£) + a) = (r(/,a)) 2 . 

Proof: The two first assertions of the lemma follow immediately from relations (16) and 
(17). The third assertion follows from the positive definitness of the matrix (D(£) + a).0 

Corollary B.2 Under the same conditions as for the previous lemma, 

i i 

(n \ 2 / „ \ 2 

C Inversion for a class of matrices 

We use the notation already introduced in Appendix B. For £ = • • • , £ n ) £ W 1 we define 
the matrix S(£) by S(£)ij = £i£ 3 for i, j = 1, • • • , n. Also we denote £~ x = (i 7 ^ 1 , • • • ,i 7 ~ 1 ),so 
that D(£- v ) = (D(i))- 1 . 
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Lemma C.l Assume £{ ^ for every i = 1, • • • , n. 

• 1/7 G R verifies 7 4' 7^ 1; ^en the matrix (D(£) — 7iS'(£)) _1 exists and verifies, 

8 = 1 

(jD( /) _ 7 5 ( i))-i = d^-i + 7 (l - 7 f>) T. 

• Reciprocally, if ^ £~i 7^ — 1; ^en £/ie matrix (D(£) + 7I)- 1 exists ans verifies, 

8 = 1 

(D(£) + 7I)- 1 = Dii- 1 ) - 7 + 7 EC 1 ) Sit" 1 )- 

Proof: Denote by R(£) the matrix R(£)ij = £j, for z,j = l,...,n. It is immediate to 
verify the identities D(£)R(£) = S(£) } R(£) 2 = (Jt ^) R{$), and = T. Hence, 

7J(£) — 7»S'(i 7 ) = D(£)(I — ~/R(£)) and therefore, it follows by direct computation (7 — 
772(f))" 1 =/| 7 (l- 7 E 4) 1 R{£). Now 

7J(f) - 7 ^(/))- 1 = TJ(f)- 1 + 7 - 7 f>) R^Di?)- 1 

which is enough to conclude. □ 
Hence, under the condition of this lemma and if 7 ^ 0, 

(D(£) - 7^(/))" 1 = D(£-^) + (7- 1 - f>) 1, 

and 

(7j(/) + 7T)- 1 = d(£-^) - [^r 1 + XX 1 ) 
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